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We investigate the conditions in which superconductivity may develop in ropes of carbon nanotubes. 
It is shown that the interaction among a large number of metallic nanotubes favors the appearance 
of a metallic phase in the ropes, intermediate between respective phases with spin-density-wave and 
superconducting correlations. These arise in samples with about 100 metallic nanotubes or more, 
where the long-range Coulomb interaction is very effectively reduced and it may be overcome by the 
attractive interaction from the exchange of optical phonons within each nanotube. We estimate that 
the probability for the tunneling of Cooper pairs between neighboring nanotubes is much higher than 
that for single electrons in a disordered rope. The effect of pair hopping is therefore what establishes 
the intertube coherence, and the tunneling amplitude of the Cooper pairs dictates the scale of the 
transition to the superconducting state. 
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I. INTRODUCTION 

During the last years there has been growing evidence 
of the existence of superconducting correlations in car- 
bon nanotubes. In the first place, the proximity effect 
has been observed in ropes of nanotubes placed between 
superconducting contacts |l],^]. Quite remarkably, it has 
been shown that the nanotubes may support supercur- 
rents below the critical temperature of the contacts [Q . 
Later on, there have been experiments directed to probe 
the superconductivity inherent to the nanotubes ||. In 
some of the ropes, a transition has been observed at a 
temperature below 1 K, with a drop of two orders of 
magnitude in the resistance down to a value consistent 
with the finite number of channels in the rope. 

From a theoretical point of view, the electronic interac- 
tions in the nanotubes have been also the subject of much 
attention . The nanotubes may show metallic prop- 
erties depending on the helicity with which the graphene 
sheet is wrapped to form the tubule. The experimental 
results H have agreed on that point with the theoretical 
predictions || . Given that the conduction takes place in 
a one-dimensional (ID) structure, it has been proposed 
that the nanotubes should be ideal systems for the obser- 
vation of the so-called Luttinger liquid behavior [p^^3[ . 
Actually, there have been measurements providing evi- 
dence of the power-law dependence of the tunneling con- 
ductance what is a signature of the mentioned 
behavior. These experiments seem therefore to probe a 
regime in which the repulsive Coulomb interaction turns 
out to be dominant in the nanotubes. 

On the other hand, the interaction with the elastic 
modes of the nanotube plays an important role in the de- 
velopment of the superconducting correlations Jlq]. The 
magnitude of the supercurrents measured in Ref. Rll is in 
some instances up to 40 times larger than expected from 
the estimate within the conventional proximity effect. It 
has been shown that the temperature dependence of the 



supercurrents is characteristic of the ID behavior of the 
system. Moreover, their values can be only explained by 
taking into account the attractive electronic interaction 
coming from the exchange of phonons, on top of the re- 
pulsive Coulomb interaction |16|. 

Recently, a microscopic model has been elaborated to 
account for the observation of superconductivity intrin- 
sic to the ropes of nanotubes [jiTj. In these systems, the 
Coulomb interaction can be significantly reduced depend- 
ing on the number of metallic nanotubes in the rope. One 
has to incorporate therefore the balance between the re- 
pulsive electron interaction and the effective attractive 
interaction coming from phonon exchange. It has been 
shown that, in the case of samples with about 300 nan- 
otubes like those displayed in Ref. j|] , the superconduct- 
ing correlations prevail in the system. The intertube co- 
herence is established mainly through the tunneling of 
Cooper pairs, which gives rise to the sup erconductivity 
in the bulk under suitable conditions |17[ ■ 

One of the aims of the present paper is to unveil the dif- 
ferent phases that arise in the competition between the 
repulsive Coulomb interaction and the attractive inter- 
action from phonon exchange in the ropes of nanotubes. 
For this purpose we will map the phase diagram of these 
systems taking the strength of the phonon couplings and 
the number of metallic nanotubes as the relevant vari- 
ables. We will show that the region where the super- 
conductivity may develop opens up for relatively large 
contents of metallic nanotubes. The phase with spin- 
density- wave correlations characteristic of a repulsive in- 
teraction is confined to the cases where there are only a 
few nanotubes in the rope. 

We assume in any event that the formation of spin or 
charge order in the bulk is prevented by the fact that 
the nanotubes have a random distribution of helicities in 
the ropes. These do not have a crystalline structure from 
the three-dimensional point of view, and the only possi- 
ble long-range order arises in the superconducting phase. 
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We will see that, between the phases with superconduct- 
ing and spin-density-wave correlations, there is a metallic 
phase with no signal of instability in any direction. This 
is the phase to be found out most likely when measuring 
ropes with a number below ~ 100 metallic nanotubes. 

We will also address the question of the maximum 
transition temperature that can be reached in the ropes. 
There have been recent experimental measures in carbon 
nanotubes of minimum diameter inserted in a zeolite ma- 
trix, from which a critical temperature of about 15 K has 
been inferred fl8|| . The small radius and high curvature 
of those nanotubes have to lead to an enhancement of the 
electron-phonon coupling. This would explain in turn the 
increase in the critical scale for superconductivity. 

One has to bear in mind, however, that the correla- 
tions in individual nanotubes provide an indirect measure 
of the superconducting state. We will show that, in the 
three-dimensional structure of the rope, the transition 
temperature is dictated by the amplitude for the Cooper 
pairs to tunnel between neighboring nanotubes. Even 
in samples with large number of metallic nanotubes, the 
opening of the superconducting phase depends at last 
on the existence of intertube coherence. The most ef- 
ficient way to increase the transition temperature may 
come actually from devising some mechanism to enhance 
the tunneling rate of the Cooper pairs between the nan- 
otubes. 

The content of this paper is organized as follows. In 
Section II we analyze the origin of the different interac- 
tions in the individual nanotubes. Section III is devoted 
to elaborate the model that takes into account the inter- 
action among the metallic nanotubes, ending up with the 
discussion of the possible phases of a rope. In Section IV 
we incorporate the effect of the tunneling of the Cooper 
pairs, performing some estimates of the transition tem- 
perature under different conditions. Finally, Section V is 
devoted to discuss some of the assumptions under which 
our general analysis applies. 



II. ELECTRON INTERACTIONS IN CARBON 
NANOTUBES 



We begin by considering the band structure of the in- 
dividual metallic nanotubes. We focus on the so-called 
armchair nanotubes, which have the structure depicted in 
Fig. " 



10 of the Appendix A. They have gapless subbands 
irrespective of the influence of the nanotube curvature 
on the orbitals. Fig. |l| represents for instance the spec- 
trum of a (10,10) armchair nanotube, which has a radius 
R I5a/ir in terms of the carbon-carbon distance a. 
We will concentrate on the case of undoped nanotubes, 
in which the Fermi level is placed at E — in the scale 
of the figure. 
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FIG. 1. Different subbands of a (10,10) armchair nanotube. 
The energy E is measured in units of the hopping parameter 
t and the momentum k is in units of the inverse of the lattice 
spacing. 

The most important feature for the study of the low- 
energy properties is the crossing of two subbands at the 
Fermi level, that we label with the indices 1 and 2 as 
shown in the figure. These subbands correspond to the 
modes with vanishing transverse momentum in the nan- 
otube, which are governed by the one-particle hamilto- 
nian 111 



U = t 



1 



1 - 2cos(V3fca/2) 

I cos (V3ka/2) 



(1) 



t being the hopping amplitude between neighboring sites. 

The amplitudes of the modes behave differently in the 
nanotube lattice depending on whether they belong to 
one or the other subband. In the case of the subband with 
bonding character, the amplitude is a smooth function 
of the position, while the modes in the other subband 
have an amplitude that alternates the sign between the 
two sublattices of the honeycomb lattice. This property 
has important consequences regarding the form of the 
electron-electron interactions as well as the coupling of 
the electrons to the elastic modes of the lattice. 

In what follows we undertake the analysis of the low- 
energy regime in which the physical properties are dom- 
inated by the branches with approximate linear disper- 
sion near the Fermi level. For this purpose, we introduce 
an energy cutoff E c below the energy of the bottom of 
the first unoccupied subband, together with its counter- 
part —E c below the Fermi level. The different branches 
are labelled by an index r = ± denoting the right or left- 
moving character, and by the Fermi point i = ± to which 
they are attached, as shown in Fig. |^. Correspondingly, 
we end up with a collection of electron fields ^ r ia {x) , to 
which one more index a will be added later on to label 
the nanotube in the rope. 
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FIG. 2. Linear branches in the low-energy spectrum of a 
carbon nanotube. 

The electron interactions that we are going to consider 
have in general the form 

fh c pk c ffcc 

dk dp I dq 

k c J —k c J —k c 

*^(p + m sj<T (p)V(k)*l ka ,(q - k)* ula ,(q) (2) 

The processes between the different branches can be cat- 
aloged according to the notation of Ref. ||. They fall 

into 12 different classes with respective couplings de- 

(i) 

noted by g\ \ the indices taking the values i = 1,2,4 
and j — 1, . . .4. The lower label pays attention to the 
Fermi points of the respective fields, following the con- 
vention that assigns i = 1 to backscattering processes, 
i = 2 to scattering between currents at different Fermi 
points, and i = 4 to scattering between currents at the 
same Fermi point. The upper label follows the same rule 
classifying the different combinations of left-movers and 
right-movers, including the possibility of having Umk- 
lapp processes (j — 3). 

We analyze next the contributions to the different scat- 
tering processes from the interactions present in the nan- 
otubes. 

Electron-electron interactions. 

We deal with the situation in which the Coulomb in- 
teraction is not screened by external gates. It is known 
that the Coulomb potential remains long-ranged in one 
spatial dimension ]20[|. The appropriate expression in 
momentum space is pl| 



V c (k) = (e 2 /4w 2 )\og\(k c + k)/k\ 



(3) 



The repulsive interaction is therefore enhanced at small 
momentum transfer, as in the processes shown in Fig. 
||. However, processes like those in Fig. EL in which 
the flavor r is changed by the interaction, are highly sup- 
pressed. This is due to the fact that, as mentioned before, 
the modes corresponding to one of the gapless subbands 
alternate the sign in the two sublattices of the nanotube. 
When the distance between the two currents that interact 
is much larger than the nanotube radius, the matrix ele- 
ment of the interaction averages to zero over the section 
of the nanotube. 
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FIG. 3. Small momentum-transfer processes which are en- 
hanced by the long-range Coulomb interaction. 
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Thus, the long-range repulsive interaction reflects in 
an enhancement of the interbranch and intrabranch cou- 

Thc flavor-changing pro- 
g^ survive in the 
form of a residual short-range interaction. Its strength, 
however, is corrected by a factor of the order of ~ 0.1 a/R 
in comparison to the nominal strength of the long-range 
Coulomb interaction |||7]]. This represents in general a 
relative reduction by two orders of magnitude, for the 
nanotubes that are typically found in the ropes. 
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FIG. 4. Flavor-changing processes with small contribution 
from the Coulomb interaction. 
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There remain the processes with momentum transfer 
of the order of ~ 2kp. These processes probe the inter- 
action at short distances. It has been shown that the 
corresponding couplings are also reduced by a factor of 
the order of ~ 0.1 a/R relative to the strength e 2 /Air 2 of 
the forward-scattering processes 00] . 

The couplings for large momentum transfer, as well as 

9% i 9^\ 9^ an( i 9^\ are marginally relevant in the 
renormalization group sense. This means that they have 
greater strength as the theory is scaled to smaller ener- 
gies, but the rate of increase starts being quadratic in 
the own couplings (|||. Thus, one has to proceed to ex- 
tremely low energies, many orders of magnitude below 
E Cl to have a significant increase of the couplings with 
respect to their nominal values || . 

Electron-phonon interactions. 

The electron-phonon interactions give rise to a re- 
tarded interaction between the electrons, which may be- 
come attractive at frequencies below the typical energy of 
the phonons. The potential for the effective interaction 
has the form 

V(k,u) = -g p , p > (k)g q>ql (k) 2 (4) 

—<JJ -+- LO k 

where Wk is the phonon energy and the g p . p i (k) are ap- 
propriate electron-phonon couplings described in the Ap- 
pendix A. The form of the interaction is dictated by 
the latter, whose dependence on the momentum trans- 
fer varies according to the kind of phonons exchanged. 

Acoustic phonons. — In this case, Uk is proportional to 
| A; | at small momentum transfer. This means that, re- 
garding the interbranch and intrabranch processes, the 
range of energies in which the effective interaction be- 
comes attractive is greatly reduced. This instance has 
been studied with detail in Ref. p2fl , looking for the 
regime in which the retarded interaction may drive the 
transition to a state with superconducting correlations. 
Taking the parameters appropriate for the nanotubes, it 
can be seen that the critical exponents estimated in this 
way are very small. Therefore, the acoustic phonons at 
small momentum transfer do not play a role in the su- 
perconductivity observed in the nanotubes. 

Moving now to the processes with large momentum 
transfer k ~ 2fcp, the phonon energy W2fc F becomes of 
the order of the Debye frequency ljd- This provides then 
a significant range in which the effective interaction be- 
comes attractive. The precise couplings to the acoustic 
phonons in the nanotubes can be found in Ref. p3| . It 
has been shown there that some restrictions arise in the 
exchange of phonons depending on the bonding or an- 
tibonding character of the subbands involved. In the 
case of longitudinal acoustic phonons, the processes can- 
not take one electron from a given subband to the other. 
Moreover, the couplings have opposite signs when the 
scattered electron is in one or the other subband. The 



exchange of phonons leads then to the effective couplings 
below the Debye frequency 

5S 1 = -ffg 1 <0 (5) 

In the case of the transverse acoustic phonons, the 
phonon-exchange processes in which the electrons remain 
in the same subband are forbidden. Processes associated 
to g^ and g[ 4 ^ in which the electron shifts from one 
subband to the other could play in principle a role, but 
there is a kinematic factor that vanishes when the incom- 
ing and outgoing modes have opposite momenta. There 
are therefore no more contributions from the acoustic 
phonons when the system is precisely at half-filling. 

Optical phonons. — The Debye frequency is of the same 
order that the energy of the acoustic phonons at 2kp, so 
that we can expect the effective couplings to be compa- 
rable to those in Eq. (||) . The different symmetry of the 
modes with bonding and antibonding character imposes 
again restrictions in the way the optical phonons can be 
exchanged, as shown in the Appendix A. In the case of 
the transverse optical phonons, the electron modes have 
to remain within the same subband, while the longitu- 
dinal phonons force the electron to shift from one of the 
subbands to the other. The effective couplings that we 
get have the following signs 

9% = ~&<0 (6) 

9% = ~9% < (7) 

& = -fl$4 < (8) 

9% = ~& < (9) 

We finally remark that the Debye frequency of the 
phonons in the nanotubes can be estimated to be between 
0.1 and 0.2 eV. This is of the same order of magnitude 
that the cutoff E c , that marks the energy below which 
the nanotubes can be viewed as simple ID systems. In 
that regime, the effective interactions coming from the 
exchange of phonons can be matched against the inter- 
action of purely electronic origin. We will adopt hence- 
forth this procedure, as long as the electronic properties 
in which we are going to be interested refer to energies 
much lower than E c . 

III. ELECTRONIC INTERACTIONS IN ROPES 
OF NANOTUBES 

We turn now to the processes that take place between 
the nanotubes in a rope. A property to be taken into 
account is that the rope is made of a disordered mixture 
of nanotubes with all kind of diameters and helicities. 
Given their random distribution, the lattices of neigh- 
boring nanotubes are not aligned in general. The conse- 
quences of this fact have been studied in Ref. p3]. As 
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long as there is a mismatch in the position of the respec- 
tive Fermi points of neighboring nanotubes, the conser- 
vation of the longitudinal momentum cannot be precisely 
enforced in the tunneling processes. At the one-particle 
level, the electronic states appear to be localized along 
the individual nanotubes, and the intertube coherence is 
strongly suppressed. 

The authors of Ref . 0] have computed the expression 
of the tunneling amplitude between nanotubes in the so- 
called compositionally disordered ropes. In the case of 
nanotubes with the same helicity and aligned lattices, 
the estimate for the tunneling amplitude is tx ~ 0.01 eV. 
When there is a mismatch dk between the Fermi points of 
neighboring nanotubes, it has been found that the ampli- 
tude for an electron hopping between them is suppressed 
by a relative factor ~ exp[— Rao(5k) 2 /A], where do is a 
parameter of the order of ~ 0.5 A. Taking a typical nan- 
otube radius R = 7 A, the average value of that factor is 
of the order of ~ 0.005. 

The presence of compositional disordered has also im- 
portant consequences regarding the kind of electronic in- 
stabilities that may develop in the rope. Depending on 
the point of the phase diagram, the individual nanotubes 
may support 2kp charge-density- wave and superconduct- 
ing correlations, or Ikp spin-density- wave correlations. It 
is clear, however, that the absence of a three-dimensional 
crystalline structure precludes the development of charge 
or spin order in the bulk of the rope. The transition to a 
superconducting state is therefore favored when the cor- 
relations are strong enough to force the symmetry break- 
ing in that direction. 

The possibility of having a superconducting transition 
in the rope of nanotubes arises from the balance of the 
effective attractive interactions due to phonon exchange 
and the repulsive interactions of purely electronic charac- 
ter. At the diagrammatic level, for instance, the correla- 
tions with s-wave order parameter are driven by the cou- 
plings g% , Qi , g% and g^\ as shown in the Appendix 
B [p5| . The obstruction to having superconducting cor- 
relations in the nanotubes comes, in principle, from the 
large contribution of the Coulomb potential to the cou- 
pling g^\ In this respect, the processes associated to g[^ 
and contribute to enhance the superconducting cor- 
relations when they are able to develop, but they play a 
secondary role in the determination of the phase diagram 
of the ropes. 

It is therefore necessary to build a model that takes 
into account the competition between attractive and re- 
pulsive interactions in the interbranch and intrabranch 
processes g^\ g^ , g± and g^\ This model can be 
solved by means of bosonization techniques. With this 
method, one boson field is introduced for each of the lin- 
ear branches below E c in each metallic nanotube. At this 
point we need to introduce a new index a labelling the 
metallic nanotubes in the rope, so that the fermion fields 



read now ^^{x). The correspondence between each of 
these and the respective boson operator is given by 

:*&V)*&(*):=P&(*) (10) 

where the colons represent normal ordering and p^ a (x) 
stands for the density operator associated to the given 
branch. 

Neglecting for the time being the tunneling between 
the nanotubes, the hamiltonian for the model, including 
a collection a — 1, .. . n of metallic nanotubes, takes the 
form 



H a = -jVF 



+ - 



di- 



ke 



—k, 



^V^(k)p^,(-k) (11) 



Jl') 



where k c is related to E c through the Fermi velocity vf, 
k c = E c /v F - 



The interaction term in (11) can be better organized 
by introducing the charge density operators 

and their spin density counterparts. We define further 
the linear combinations 



p%(x) = ±(p%\(x)+p%l(xj) (12) 



(13) 

p ( 2 f !(x)=P i : ) - P (x)+p { X(x) (w) 



» 



Obviously, p\'{x) is made of the modes with bonding 

(a) 
2p 



character, while p2p( x ) contains those with antibonding 



character. 

We recall that, in the case of the attractive interaction 
from the exchange of phonons, the sign of the coupling 
depends on the bonding or antibonding character of the 
currents involved. Thus, a great simplification is achieved 
by introducing the operators 



pf p {x) = ±{fi p \x)+pf p \x) 



(15) 
(16) 



The Coulomb interaction acts on the symmetric combi- 
nations p+p(x), while the effective attractive interaction 

only couples the antisymmetric combinations p_ p {x). 
The hamiltonian ( |ll"l) can be rewritten in the form 

#o = \vf J ' dk £ : pH(k)pH(-k) : 



b 



+gY,~P { -l(k) (&(-k) 



(17) 
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where we have defined 
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(2) (4) 
92 =94, = 



(4) 

-92 = 



(2) 

-9l 



(18) 



After the decoupling of the Coulomb and the effective 
attractive interaction, we observe that the latter affects 
to a number n of the An density operators, while the 
Coulomb interaction acts only on the channel of the total 
charge density. That is, the hamiltonian (|l7]) can be com- 
pletely diagonalized by changing variables to the sym- 
metric combination of the charge density operators of all 
the metallic nanotubes. This is a consequence of the long- 
range character of the Coulomb interaction, that allows 
us to assume that, in the interbranch and intrabranch 
processes, the different nanotubes interact between them 
with the same strength. 

The balance between the long-range repulsive interac- 
tion and the effective attractive interaction is possible as 
the former turns out to be greatly reduced for large values 
of n. The attraction from phonon exchange is otherwise 
an intratube effect that does not depend on the number 
of nanotubes. It is then interesting to draw the phase 
diagram of the ropes in terms of the number of metallic 
nanotubes and the strength of the attractive interaction. 
At large values of n and sufficiently large values of \g\, 
there must exist a phase where the attractive interaction 
prevails, with the development of superconducting cor- 
relations in the individual nanotubes. For low values of 
n and small |g|, the electronic properties have to be dic- 
tated by the Coulomb interaction, with the appearance 
of the spin-density- wave correlations characteristic of the 
Luttinger liquid behavior. 

The different phases can be identified by looking at the 
correlators of the model. The bosonization techniques 
allow to compute them by using the correspondence be- 
tween the fermion fields and the respective boson opera- 
tors @0 



*M a (x) = e X p(-i<!>™ tr (x) 



(19) 
(20) 



where d x ^ a (x) = 2np^ a (x) . 

We deal first with the propagator {x, t) of the 
Cooper pairs in the individual nanotubes. This object 
factorizes into the different channels that appear after 
diagonalizing the hamiltonian 

Dj$(x,t) = K°it(M)*--J (*, ^1,(0,0)^(0,0)) 



3n-l 



C sc (x,t)Y[N sc (x,t) Yl F(x,t) 



(21) 



The first factor corresponds to the channel of the to- 
tal charge density of the metallic nanotubes, while N sc 
stands for the contribution of the antisymmetric combi- 
nation of the charge in the two subbands. The rest of the 
factors corresponds to the channels with no interaction. 



The different factors in Eq. ( |21| ) can be computed in 
terms of the respective boson propagators, what leads to 
expressions of the form 



X sc (x,t) = exp ( -— / dk 
\ 2n J 



/i(fc) k 



(1 - cos(kx) cos(v F kt)) (22) 



In the case X sc (x,t) = C sc (x,i), we have /i(fc) = 
+ 8nV c {k)/v F and v F {k) = v F /u(k). For 
X sc (x,t) = N sc (x,t), we have the constant parameters 
u = 1/ ^/l — A\g\/Trv F and v F = v F j \x. 

The same technique can be applied to the computation 
of the propagator D^ w (x 1 t) of the spin-density-waves 
along the nanotubes. The correlator is given by 



^dt(M) = (*^(x,0*^;(x J *)*L 8 i\(0 > 0)*V+T( ' )) 

n 3n— 1 

= C S t)l[N sdw (x,t) Yl F(x,t) (23) 
1 1 

Each of the factors in Eq. ( p3|) has now the representation 

X sd w(x,t) = exp [ / dk 
\ 2n J 



( a )t In n\jiW 



(1 - cos(fcr) cos(v F kt))j (24) 

with the same correspondence between fi(k) and C, N 
and F as for the propagator of the Cooper pairs. 

The correlators in the individual nanotubes cannot dis- 
play an instability for any finite value of the frequency 
p3| . The way to discern the tendency to the formation of 
long-range order is to determine whether the correlations 
are enhanced at large distances over the values in the ab- 
sence of interaction. When this happens, the Fourier 
transform of the corresponding propagator at zero tem- 
perature diverges at zero frequency and momentum. At 
nonvanishing temperature, the propagators remain finite 
but the intertube coupling may give rise to the break- 
down of the symmetry if the correlations in the nanotubes 
are sufficiently enhanced. We will address this question 
in the next section. The computation of the propagators 
( pl| ) and (|H|) can be extended to the case of temperature 
T^O, just by inserting the factor 1+2 / \exp(v F \kUT) — 1] 
in the integrand of expressions such as (|22| ) and (|24|) [0 . 

At this point, we stick to the model at zero tempera- 
ture and map the regions with enhanced superconduct- 
ing or spin-density- wave correlations. The propagators 
( pl| ) and ( p3| ) do not show a perfect scaling behavior at 
large distances, due to the momentum dependence of the 
Coulomb potential. We may take however an effective 
value fco of the momentum in the infrared to approxi- 
mate the behavior of each propagator by a power-law 
with a constant exponent. 
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In the case of D 
takes the form 



(0) 



(x,t), the decay at large distances 



D£\x,0)~l/x^ 
where the anomalous scaling dimension is given by 



(25) 



2n/i(fco) 2fi 



1 

2n 



Similarly, (x, t) has a large-distance behavior 



D 



(0) 



'sdw(^'°) ~ l l X 

with an anomalous exponent 
H(k 



2+8 



8 



2)i 



2 



1 

2n 



(26) 



(27) 



(28) 



We have represented in Fig. |5| the different phases that 
arise by varying the number n of metallic nanotubes and 
the strength \g\ of the intratube attractive interaction. 
For the Coulomb interaction, we have taken the values 
2e 2 /n 2 VF — 1.0 and ko — 10~ 3 fc c , which are appropri- 
ate for typical experimental samples Jl(|. At sufficiently 
large values of n, a phase with superconducting correla- 
tions opens up, characterized by negative values of the 
anomalous dimension 7. This phase is placed above the 
upper full line in the diagram. Another phase with spin- 
density- wave correlations shows up, below the lower full 
line in the diagram. This region corresponds to the points 
characterized by negative values of 5. 
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FIG. 5. Phase diagram showing the regions where the su- 
perconducting (SC) correlations and the spin-density-wave 
(SDW) correlations prevail, in terms of the strength \g\ of the 
effective attractive interaction and the number n of metallic 
nanotubes. 

We observe also the existence of an intermediate re- 
gion in which neither 7 nor 5 are negative. This rep- 
resents a phase with intrinsic metallic properties, where 
no correlations are enhanced and no order would develop 
irrespective of the value of the intertube coupling in the 



rope [g6[. In our model, this phase arises from the bal- 
ance between the Coulomb and the effective attractive 
interaction. The former remains long-ranged, but the in- 
teraction among a large number of metallic nanotubes is 
what favors the appearance of the metallic phase, as a 
bridge to the development of the superconducting corre- 
lations. 



IV. SUPERCONDUCTING TRANSITION IN 
ROPES OF NANOTUBES 

In the preceding section we have neglected the effect 
of tunneling between metallic nanotubes, relying on the 
smallness of the amplitude for that process. We have 
seen however that, in the absence of intertube hopping, 
the model cannot develop any instability at finite tem- 
perature. In this respect, the upper full line in Fig. ||has 
to be considered as a first approximation to the bound- 
ary where the superconducting phase opens up in real 
samples. The tunneling amplitude between neighboring 
nanotubes in a rope is actually what dictates the tran- 
sition temperature to the superconducting state, as we 
study in what follows. 

As mentioned before, the single-particle hopping be- 
tween neighboring nanotubes is strongly suppressed in a 
disordered rope. The tunneling amplitude between nan- 
otubes with different helicities can be estimated to be of 
the order of ~ 0.5 x 10~ 4 eV. Its smallness stems from 
the obstruction to having precise momentum conserva- 
tion in the hopping between the misaligned lattices of 
the nanotubes 0|. In these conditions, the tunneling 
of Cooper pairs turns out to be a much more important 
effect. 

The probabilities of the single-particle hopping and the 
pair hopping can be compared by measuring the respec- 
tive amplitudes with respect to the energy cutoff E c . 
We can take this quantity of the order of ~ 0.1 eV, 
which is the scale below which the nanotubes in a rope 
are seen as ID objects. Thus, we find that the single- 
particle hopping has a relative weight of the order of 
~ 0.5 x 10~ 3 . The amplitude for the tunneling of Cooper 
pairs can be computed roughly as the square of the am- 
plitude frr ~ 0.01 eV, so that the relative weight A2 of 
that process becomes of the order of ~ 10~ 2 . 

The model of the preceding section has to be corrected 
then by adding to the hamiltonian (0) the term describ- 
ing the tunneling of the Cooper pairs 



ck c rk c 

,/, 1 dk I dp 



fee 



dp' 



H2 — ^ (A2)ofc 

(a, 6) J -k c J-k c J-k c 

*lf(k +p)*L a ^ ii (-p)*L b L jU (-p')*S ) T (A : +p>) (29) 

where the sum runs over all pairs (a, b) of nearest- 
neighbor metallic nanotubes. 
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Given that the model without pair hopping is exactly 
solvable, we can compute the scaling dimension of the 
operator in ff 2 . We can check in this way whether this 
is a relevant perturbation or otherwise it fades away at 
low energies. The naive dimension of the four-fermion 
operator is 2, and the deviation from that value can be 
obtained from the large-distance behavior of the correla- 
tor in Eq. (pl|). The anomalous scaling dimension coin- 
cides with the exponent given in Eq. (|2^). Therefore, we 
observe that it is negative in the phase with supercon- 
ducting correlations, what points at the enhancement of 
the pair-hopping perturbation at low energies. 

The values that can be obtained for 7 are rather small 
anyhow. Taking again an effective value of fco = 10 _3 fc c 
, we find 7 w —0.02 for n — 60 and 4\g\/irvp = 0.6, and 
7 w -0.16 for n = 120 and A\g\/irv F = 0.8. In the lat- 
ter case, for instance, the anomalous scaling leads to an 
enhancement by a factor that, four orders of magnitude 
below E c , can be at most of the order of (10 -4 ) 7 ~ 4. 
Such a scaling does not affect our estimate of the order 
of magnitude for the relative weight A 2 . 

As long as we have the small parameter A2 in the sys- 
tem, we may include the effects of pair hopping by sum- 
ming up all the processes in which a Cooper pair, propa- 
gating along a nanotube with the amplitude ( pl| ) , tunnels 
to a neighboring metallic nanotube. Let us denote by l a 
the position of nanotube a in the transverse section of 
the rope. We will assume that the metallic nanotubes 
are dense in the collection of nanotubes of the rope |27| . 
Then, we may write the propagator of the Cooper pairs 
between the metallic nanotube a and the metallic nan- 
otube b as a function D(l a , 4; x, t). This object is related 
to Dsc\x, t) through the self-consistent equation repre- 
sented in Fig. ||, in the approach that takes H2 as a 
perturbation to the hamiltonian (p"7|). 



D 
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FIG. 6. Self-consistent diagrammatic equation for the 
propagator D of Cooper pairs along the rope. 

For the sake of simplifying the calculation, we will sup- 
pose that the positions l a form a periodic arrangement in 
the transverse section of the rope. Then we can take the 
Fourier transform of these variables as well as of the dis- 
tance x along the nanotube. The equation for the Fourier 
transformed propagator D(q;k,u>k) reads 

5(g; k, Lu k ) = (fc, LO k ) + D<$ (fc, w fc ) A 2 (q)D(q; k,u k ) 

(30) 

The measure of the condensation of Cooper pairs in 
the rope is given by the propagator at zero frequency and 



momentum D(0; 0, 0). This accounts for the propagation 
of a Cooper pair from a metallic nanotube to the rest in 
the rope. From Eq. (BCf) we obtain the expression 



0(0; 0,0) = 



^ c 0) (o,o) 



1-A 2 (0)CL 0) (0,0) 



(31) 



According to the conventional interpretation, the tran- 
sition to the superconducting state is given by the point 
at which D(0; 0, 0) develops a pole. This may happen 
only when Dsc^(0,0) becomes large enough at low tem- 
peratures. In the superconducting region of the phase 
diagram, the correlations grow large in the limit of van- 
ishing temperature, as shown in Fig. |?]. The only limi- 
tation to develop a real divergence is placed by the finite 
length of the system, as we discuss later on. 
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FIG. 7. Plot of the propagator at zero frequency and 
momentum versus T/E c , for 2e 2 /(-K 2 vp) = 1.0 . The dashed 
line corresponds to the case n = 1 and g = 0, and the solid 
lines to n = 100 and the respective values (from top to bot- 
tom) 4\g\/(nv F ) = 0.75,0.5,0 . 

The parameter that plays the major role in setting the 
values of the transition temperature is the weight A 2 (0) 
for pair hopping at zero transverse momentum. One can 
look, for instance, for the points of the phase diagram 
with transition temperature T c = 10 _3 i? c . This value 
corresponds to a temperature of w 1 K for E c = 0.1 eV. 
Taking the intermediate value A 2 (0) = 0.025, the points 
form the boundary represented by the dashed line in Fig. 
|. The curve has the same shape that the boundary of the 
superconducting phase determined from the expression 
( p6| ) of the anomalous dimension. However, we see that 
the region with T c > 10~ 3 E C above the dashed line is 
sensibly smaller than that of the whole superconducting 
phase. 

We have represented in Fig. || the contour lines for 
different critical temperatures in the space of the pair 
hopping parameter A 2 (0) and the number n of metallic 
nanotubes, fixing the coupling of the attractive interac- 
tion at 4\g\/-KVF — 0.75. We observe that a slight change 
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in the value of the tunneling amplitude leads to a signifi- 
cant increase in the transition temperature. Considering 
ropes with larger content of metallic nanotubes may also 
help to enhance T c , although the figure shows that very 
high values of n have to be reached to find a sensible 
variation. 

0.06 t 1 



0.04 



A; 



0.02 - 



0.00 
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10 E c 



0.5x10 E„ 



100 



10 E c 
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The constraint on the superconducting correlations 
due to the finite-size scaling may have been observed in 
the experiments reported in Ref. ||. It is shown there 
that the large drop in the resistance is present in two 
samples with respective lengths of 1.6 /im and 1 /im. 
The effect is absent in a third sample which is 0.3 /im 
long. Assuming that the typical transition temperature 
for these samples is around T c ~ 0.5 K, the relative small 
length of the third sample would explain the absence of 
superconductivity. The minimum transition temperature 
that could be supported in that case is of the order of 
T c ~ 0.1 v F /L ~ 10~ 4 eV - 1 K. The same argument 
provides a definite check of the model elaborated in this 
paper since, according to it, a transition with critical tem- 
perature T c = 0.1 K, for instance, should not be present 
in the samples below 1 /im. In general, it has to be true 
that any sample with length L cannot have a transition 
at temperatures lower than 0.1vf/L, as it is observed in 
the samples considered in Ref. 0] . 



FIG. 8. Contour lines for the critical temperature T c de- 
pending on the pair-hopping parameter A2 (0) and the number 
n of metallic nanotubes. 

Finally, we remark that the finite length of the nan- 
otubes imposes a limit on the strength of the correla- 
tions. The existence of a length scale L spoils the scaling 
of the model and, therefore, the approximate power-law 
behavior of the correlators. Thus, the divergence of the 
propagators at zero frequency and momentum is cut off 
in practice at a temperature scale which is between vp/L 
and one order of magnitude below that value. This ef- 
fect is illustrated in Fig. ^, which displays the behavior 
of 5L 0) (0,0) for ropes with different number of metal- 
lic nanotubes n = 40,100,400,1000 and finite length 
L = 1000/fc c . 
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FIG. 9. Logarithmic plot of at zero frequency 

and momentum versus T/E c , for 2e 2 / \-k 2 vf) = 1.0 and 
4\g\/(irvF) = 0.75. From top to bottom, the curves 
correspond to different numbers of metallic nanotubes 
n = 1000, 400, 100 and 40 . 



V. DISCUSSION 

In this paper we have shown that superconductivity is 
a plausible effect in the ropes of nanotubes. The ropes 
with greater number of metallic nanotubes have in gen- 
eral larger superconducting correlations. We have seen 
that the long-range Coulomb interaction is reduced very 
effectively in ropes with 100 metallic nanotubes or more. 
On the other hand, the coupling to the elastic modes 
within each nanotube provides the attractive interaction 
leading to the electronic pairing. We have shown that the 
low-momentum optical phonons have suitable properties 
to balance the effect of the repulsive Coulomb interac- 
tion. 

We have dealt with a model providing an exact de- 
scription of the competition between the Coulomb in- 
teraction and the effective attractive interaction in in- 
terbranch and intrabranch processes. Proceeding in this 
way, we have disregarded the effect of other processes 
coming from phonon exchange. These contribute to the 

backscattering and Umklapp couplings g± , 9^, 9% and 

(3) 

g 2 ■ Only the first and the third kind of processes play a 
role in the development of the superconducting correla- 
tions, as analyzed in the Appendix B. These coupling are 
marginally relevant in the renormalization group sense 
but, as we have already remarked in the paper, the ef- 
fects derived from renormalization group scaling turn out 
to be very soft down to the energies where to look for a 
superconducting transition. 

The significance of the backscattering processes is 
found in that they determine the symmetry of the order 
parameter whenever the system becomes superconduct- 
ing. According to the analysis in Section II, both 
and 02 correspond to an effective attractive interaction. 
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Then, from inspection of the different order parameters 
considered in the Appendix B, we conclude that singlet 
pairing is enhanced by the backscattering interactions. 
We observe also that the s-wave symmetry with positive 
amplitude in all the branches is favored over more ex- 
otic possibilities like the d-wave symmetry of the order 
parameter. 

In our description of the mechanism of superconduc- 
tivity, we have taken into account the fact that the ropes 
are made of nanotubes with a random distribution of he- 
licities. This prevents the development of any ordered 
charge or spin structure from the coupling of the nan- 
otubes. Moreover, this kind of disorder implies a strong 
suppression of the single-particle hopping and that the 
relevant tunneling process is given by the hopping of 
Cooper pairs between neighboring nanotubes. The am- 
plitude for that process is in general very small, which 
explains the relatively low transition temperatures mea- 
sured experimentally. 

In any event, the effect of pair hopping has the most 
direct influence on the transition to the superconducting 
state. The disorder present in the samples used in the 
experiments has to be very high, and its effect should be 
studied on more quantitative grounds to have an estimate 
of the maximum transition temperature reachable in the 
ropes. In this sense, good prospects should exist to in- 
crease the transition temperatures, either by suitable re- 
finement of the experimental samples or by intercalation 
or modification of the internal structure of the ropes. 

Fruitful discussions with F. Guinea and A. Kasumov 
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supported by CICyT (Spain) and CAM (Madrid, Spain) 
through grants PB96/0875 and 07N/0045/98. 



APPENDIX A: SELECTION RULES FOR THE 
COUPLING TO OPTICAL PHONONS 

The different symmetry of the modes with bonding and 
antibonding character leads to definite relations between 
the various electron-phonon couplings. This can be ob- 
served in the results of Ref . (23) , where the full expres- 
sions of the couplings to the acoustic phonons in armchair 
nanotubes have been obtained. In this appendix we ex- 
ploit the mentioned symmetry to get the corresponding 
relations in the case of the optical phonons. 

The use of the tight-binding approximation is appro- 
priate for the carbon nanotubes p^] . The electron- 
phonon coupling can be written in terms of the polariza- 
tion vector e s (k) depending on site s and the amplitudes 

of the incoming and outgoing modes uf^ (k) and uf, ' (k 1 ) 
in the respective subbands p and p', taking the form 

9 P , P '(k,k') = ^ T7 J ]T U W*(k)t#V) 



(e s (k - k') - e s ,(k - k'))-VJ(s, s') (Al) 

In the above expression, the sum is restricted to nearest- 
neighbors of the atoms in the unit cell of the nanotube, 
J(s, s') is the matrix element of the atomic potential con- 
necting orbitals at s and s', tOk is the phonon frequency 
and /z is the mass per unit length. 

Let us consider first the case of optical phonons with 
longitudinal polarization. The polarization vector de- 
pends only on the distance z measured along the nan- 
otube, but it has opposite amplitudes in the two sublat- 
tices shown in Fig. [l(] so that 

e s (k) = z exp(ikz s ) for black points (A2) 
e s (k) — — z exp(ifcz s ) for white points (A3) 

i being the unit vector in the axis direction. 




Z 



FIG. 10. Scheme of the unit cell of an armchair nanotube, 
represented between the vertical dashed lines. The arrows 
correspond to the displacements for a longitudinal optical 
phonon. 

It is not difficult to see that, if the modes belong to 
the same subband p = p', the different terms in the sum 
of the expression ( |Al| ) cancel among themselves. Then, 
we have that in the case of longitudinal optical phonons 

9i,i(k,k')=g 2t2 {k,k') = (A4) 

When the incoming and outgoing electron modes are in 
different subbands, the electron-phonon coupling does 
not vanish. Using the fact that the amplitude of the 
modes with antibonding character changes its sign under 
the exchange of the two sublattices, we obtain the result 

gi,2(k,k') = -g 2 ,i(k,k') (A5) 

Moving now to the case of transverse optical phonons, 
we have a picture like that shown in Fig. |ll|. The polar- 
ization vector is 

e s (k) = 8 cxp(ikz s ) for black points (A6) 
e s (k) = —9 evp(ikz s ) for white points (A7) 
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9 being a unit vector tangential to the nanotube. When 
the incoming and outgoing electron modes belong to the 
same subband, we have now 

g hl (k,k , ) = -g 2 , 2 (k,k') (A8) 

If the modes are in different subbands, the terms in the 



expression (Al) cancel out in pairs and we have 

Si |2 (fc J /E / )=52 J i(fc,fc') = (A9) 



The symmetry of the superconducting correlations can 
be obtained by using a diagrammatic analysis of the 
different response functions. In the case of singlet and 
triplet pairing, the contributions to the correlator have 
the structure depicted in Fig. [l2L with respective up- 
per and lower signs. Then, the enhancement of the sin- 
glet pairing response function is driven by the combina- 
tion of renormalized couplings g% + g^ + + 9^\ 

while in the case of triplet pairing the combination is 

12) (i) , (l) (2) 



9i -9\ 



92 



9i 



FIG. 11. Same scheme as in Fig. [j^, but representing the 
displacements for a transverse optical phonon. 



We observe that the situation is reversed with respect 
to the case of the acoustic phonons, in which the relations 
(A8) and (|A9D a Pply to the longitudinal branch while 
(A4) and (|A5|) hold for the transverse branch. 



APPENDIX B: PAIRING SYMMETRY OF THE 
SUPERCONDUCTING CORRELATIONS 

There are several response functions that give a mea- 
sure of the superconducting correlations in the nan- 
otubes, each of them corresponding to a different sym- 
metry of the pair wavefunction. They have the general 
form 



R{x,t) = (0(x, t)0t(o,O)) 



(Bl) 



In the case of singlet and triplet pairing, the operator 
0(x, t) corresponds respectively to the upper — sign and 
the lower + of the expression 



0(x, t) = * ++T *__j t- t 
+ v f' — hT^H — 1 "F ^ — hl^H — T 



(B2) 



There are also more exotic possibilities like the d-wave 
symmetry of the Cooper pairs, corresponding to 



-(*-+?*+-;- *-+;*+-t) 



(B3) 




+ 



+ 





+ 




FIG. 12. Contributions to the response functions for sin- 
glet and triplet pairing, where the filled circles represent the 
different renormalized vertices. 

For the response function with <i-wave symmetry of 
the Cooper pairs, the structure of the contributions is 
represented in Fig. |l3|. The enhancement is controlled 
now by the combination of couplings g^ + g^ 1 




++t 



+ 



— \ 






FIG. 13. Contributions to the response function for d-wave 
pairing, with the same representation as in Fig. |l| 
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